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Abstract 

The six Painleve transcendants which originally appeared in the studies of ordinary 
differential equations have been found numerous applications in physical problems. The 
well-known examples among which include symmetry reduction of the Ernst equation 
which arises from stationary axial symmetric Einstein manifold and the spin-spin cor- 
relation functions of the two-dimensional Ising model in the work of McCoy, Tracy and 
Wu. 

The problem we study in this paper originates from random matrix theory, namely, 
the smallest eigenvalues distribution of the finite n Jacobi unitary ensembles which 
was first investigated by Tracy and Widom. This is equivalent to the computation of 
the probability that the spectrum is free of eigenvalues on the interval [0, t] . Such 
ensembles also appears in multivariate statistics known as the double- Wishart distri- 
bution. 

We consider a more general model where the Jacobi weight is perturbed by a discon- 
tinuous factor and study the associated finite Hankel determinant. It is shown that the 
logarithmic derivative of Hankel determinant satisfies a particular a— form of Painleve 
VI, which holds for the gap probability as well. We also compute exactly the leading 
term of the gap probability as t — > 1~ . 



1 



1 Introduction and statement of results 



The discovery by Tracy and Widom [29], [30] in the 1990's of probability laws that govern 
the extreme eigenvalues of certain Hermitian random matrices triggered a renewed interest 
in random matrix theory and has lead to an explosion of activity in recent years; cf. [3] for a 
recent account regarding the application of random matrix theory to wireless communications 
and multi-variate statistics. In particular, Painleve transcendants play an important role in 
these studies. 

Painleve transcendants refer to solutions of six nonlinear second order ordinary differential 
equations of the form 

y" = R(y,y',t). 

Here ' denotes -| and R is a function rational in y, y' and analytic in t. All Painleve tran- 
scendants enjoy the celebrated Painleve property which asserts the only movable singularities 
are poles. 

There are broadly speaking two well-known physical applications of the Painleve equa- 
tions. One is the reduction of the integrable Ernst- Weyl equations to Painleve VI by Cosgrove 
[32] and the reduction of the integrable Loewner-Konopelchenko- Rogers system [T7] to a ma- 
trix form of Painleve VI by Schief [27J . For a review on the reduction of the Ernst equation, 
we refer to [21]. The other application is to the two-dimensional Ising model where McCoy, 
Tracy and Wu [22] showed that the spin-spin correlation function expressed as the Fredholm 
determinant of certain integral operator is a particular Painleve III. 

In this paper, we are also concerned with Fredholm determinant of a certain finite rank 
integral operator, which is equivalent to the determinant of a Hankel matrix generated by a 
weight of the form 

x a (l -xf(A + B9(x-t)), are [0,1] (1.1) 

where t G [0, 1] , a > , (3 > and #(•) denotes the Heaviside step function. 

We take the shifted Jacobi weight, x a (l — x) , to be the "reference" weight Wq and call 
w wj a perturbed Jacobi weight, where wj(x;t) = A + B9(x — t) is the "jump" function. 
The constants A and B are chosen in such a manner that A > and A + B > , which 
ensures the real Hankel form 

n-l 

^ CjCj llj+jjt), CjfER, 2 = 0,1,..., 
i,j=0 

is positive definite and the Hankel matrix defined below 

Tin := (Vi+j(t))™ j=0 , 

is invertible, where 

Vi{t) := / x a+i {l-xf[A + Be{x-t)]dx, (1.2) 
Jo 
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for i = 0, 1, 2, ... . An evaluation of the integrals in (jl.2p gives 

B 



fii(t) = AT(i + a + 1) + 



1 _ t f+P f+a 2 Fx(-a - i, 1; 2 + /?; 1 - 1/t), 



where T(x) denotes the Gamma function and 2 -Pi is the hypergeometric function; see [T]. 
It is a well-known fact that the Hankel determinant 



vn-l 



D n (t) = detft n = det (^+i(*))" i=0 
has the following multiple integral representation [28] 



[0,1] 



[A n (s)] 2 Y[ x% (1 - ac fc ) t)dx k , 



;i.3) 



:i-4) 



fc=i 



where 



AJx) 



n 1 

l<i<j'<n 



Hence, we can also view D n (t) as the partition function for the random matrix ensemble 
with the joint eigenvalue distribution 



which is of interest in statistical physics. 
Let 



\\xf{l- Xl f (A + Be( Xl -t))d Xl , 



d 



H n (t) :=t{t-l)-\nD n (t). 



1.5) 



It is the aim of this paper to produce a second order differential equation (non-linear) that 
is satisfied by the above logarithmic derivative of D n (t) . Our main result is the following: 



Theorem 1 The a function, defined as a translation in t of H n , 

a{t) := H n {t) + dit + d 2 

with 

d 1 = - n (n + a + /3)- - — 
1 



;i.6) 



do 



2n(n + a + (3) + j3(a + f3) 
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satisfies the following Jimbo-Miwa-Okamoto flH$ , WM) o~—form of Painleve VI 
a'{t{t - l)a"\ + 1 2a' {ta' -a)- a' 2 - z/iZ/ 2 z/ 3 z/ 4 



= (a> + v 2 )(a> + v 2 2 )(a> + v 2 3 )(a> + v 2 4 ) (1.7) 

with 

a + (3 P — a 2n + a + ft 

and the initial conditions 



vi = — — , v-i = — ~ — , ^3 = ^4 



cr(0) = d 2 , o-'(O) = di. 

Note that the constants A and B do not appear in the differential equation (11.71) . 

In the situation where A = and B — 1 , the quotient D n (t) / D n (0) stands for the gap 
probability that the interval [0, t] is free of eigenvalues or, equivalently, the least eigenvalue 
of (the finite n) Jacobi unitary ensembles is > t. By a fundamental result of Gaudin and 
Mehta |23j, the gap probability can be formulated as the following Fredholm determinant: 



ggft) _ dct(I K W) 
D n (0) ~ 1 



The integral operator acts on functions / G L 2 (0,t) as 

Jo x y 

where constant explicitly known and 4> n {x) is a constant multiple of x Q//2 (1 — 

with Pn"'^ (x) being the Jacobi polynomials. Our next result deals with the 
asymptotic behavior of gap probability as t — > 1~ : 

Theorem 2 Suppose that A = and B = 1 , we have 

Dn(t) 2 2na _ ,N„(„ +/ ?) Q g x 

D n (0) F(a,(3)K(a,f3,n) { ' 1 ' J 

as t — > 1~ , where 

F( _ r((a + g + l)/2)G 2 ((a + (3 + 1)/2)G 2 (1 + (a + (3)/2) 

{,P) G{a + (3 + l)G{a + 1)G{(3 + 1) ' { ' 

, , _ G{n + l)G(w + a + l)G(n + + l)G(w + « + /? + !) 

(«, (J, n) - G2{ ^ + (a + ^ + 1)/2)G2(n + ! + (a + /3)/2)r(n +(« + /?+ l)/2) ' ( ' j 

ifere, G(x) denotes the Barnes G- function: 

oo 

G(x + 1) = ( 27r )*/V^ +1 )+^/ 2 J] [(l + ^V*+- 2 /(2n) 

n=l 

where 7 t/ie Euler-Mascheroni constant. 
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We mention here that the so-called Dyson- Widom constant; the coefficient of (1 — t) n( - n+ ^ 
of (1.8) was conjectured in [TB] for the gap problem, but with the more general hypergeo- 
metric kernel. However that result cannot be applied here since the parameters there and 
those of our problem are different. 

The differential equation satisfied by the logarithmic derivative of D n (t) when A = and 
B = 1 is one of the problems first investigated by Tracy and Widom in [31] using operator 
approach of [29J and [30]. This method brings in the resolvent of K^l and distributional 
kernel, gives rise to a third order ordinary differential equation. It was not clear how a first 
integral could be found, as one expects to find a second order equation. Recently, a first 
integral was found in [18] for the Gap problem for hyper geometric kernel. 

Later, Haine and Semengue [15], made use of the "infinitely-many-times" approach of 
Adler, Shiota and van Moerbeke [2], popular amongst the experts in integrable systems, to 
investigate the gap problem in the Jacobi polynomials ensemble. This approach involves the 
notion of vertex operators, an important tool in Sato's theory. Another third order ordinary 
differential equation was found. Eliminating the third derivatives from both equations, 
Haine and Semengue, found that the resulting second order ordinary differential equation is 
a particular Painleve VI. 

The appearance of Painleve VI could perhaps be guessed at based on the following con- 
siderations: In the original or "unperturbed" problem, where B = 0, the Jacobi polynomials 
satisfies a second order ordinary differential equation with regular singular points at 0, 1, 
and oo. It is to be expected that "switching" on the discontinuity (B > ), an extra regular 
singular point is introduced at t . Because the resulting ordinary differential equation has 
four regular singular points, one of which moving, according to the isomonodromy theory of 
Jimbo and Miwa [IB], one expects to find Painleve VI and allied functions appearing in the 
coefficients of the ordinary differential equation. However, it is not immediately clear how 
the correct parameters may be identified, and it is also not clear which quantities related to 
the orthogonal polynomials is that particular Painleve VI. 

Because of the discontinuity of the weight at t and the t dependence of the moment 
sequence {/ij(t)}, i = 1,2,..., the coefficients of polynomials {P n (z)} orthogonal with 
respect to w wj also depend on t. We normalize our monic P n (z) to be 

P n (z) =z n + p x {n,t)z n - 1 + ... + P n (0). 

It turns out that 

HJt) =t(t- 1) — In DJt) = -(2n + a + /3)p 1 (n,t) - n(n + a); 
at 

see Theorem in Section [31 This, together with (11.61) . implies that our a— function is up to a 
translation in t essentially p 1 (n, t) , the coefficient of z n ~ l of monic orthogonal polynomials 
P n (z) • Comparing with the Tracy- Widom approach, we see that their resolvent evaluated 
at the same point is also interpreted here as p x (n, t) . 

The approach in this paper is based on the ladder operators of orthogonal polynomials 
and the associated compatibility conditions (Si), (S2) and (S' 2 ) , the detailed derivation 
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of which can be found in [I]. We find through the ladder operators there are four auxiliary 
quantities x n (t), y n (t), R n (f) and r n (t). This appears to create complications. But the 
compatibility conditions, as we shall see later, ultimately provides a determination of these 
through a system of difference equations. 

For a comparison between the ladder operator theory and the isomonodromy theory of 
Jimbo and Miwa [16] in a particular Hermitian random ensembles, see [10] . There the poly- 
nomials are orthogonal with respect to a "singularly" deformed Laguerre weight; x a e~ x e~ s ^ x , 
over M + , and one finds a Painleve III. The reader is also referred to [2E] for a comparison 
between the Tracy- Widom theory and that of Adler, Shiota and van Moerebeke. 

In the development of what follows, we would like to remind the reader to keep in mind 
that the recurrence coefficients a n {t) and (3 n {t) (see (12.31) below) play a central role and 
attention should be focused on expressing these in terms of the auxiliary quantities mentioned 
above. The reader will see that the auxiliary quantities appear naturally in the ladder 
operator approach. We will show that one of the auxiliary quantities x n (t) , intimately 
related to the recurrence coefficients, is up to a linear fractional transformation a particular 
Painleve VI. 

This rest of the paper is organized as follows. In the next section, we introduce the 
notations to be used throughout the paper and derive two Toda-type equations; these are 
differential- difference equations involving a n {t) and (3 n {t) . In Section [31 a system of non- 
linear difference equations involving the auxiliary quantities are found and are used to obtain 
expressions of the recurrence coefficients amongst other things. Based on these facts, we 
prove Theorem [1] and show that x n (t) -one of the auxiliary quantities-is actually a particular 
Painleve VI in Section HI We conclude this paper in Section with the proof of Theorem [21 

2 Preliminaries: notations and time evolution 

The purpose of this section is to derive two universal Toda-type equations which involve 
the recurrence coefficients and the auxiliary quantities R n {t) and r n {t) . We keep w quite 
general. 

The general To da-hierarchy, obtained through the multiplication of wo by exp ^— Y^JLq tj x ^ 
can be found for example, in [21], [2] and [Hj . See also [15] for a discussion of the approach 
of [2] in relation to Sato's theory. 

To begin we consider polynomials {Pi(x)} orthogonal with respect to the weight 

w (x)(A + B6(x-t)) 
on [0,1]. The weight wq will be known as the "reference" weight, and 

wj(x;t) := A + B6(x - t) 

as the "jump" factor. 
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The orthogonality condition is 

P i (x)P j (x)w (x)w^x;t)dx = hi(t)5ij, fn(t) > (2.1) 



Jo 



for i — 0, 1, ... , and the t dependence through wj, induces t dependence on the coefficients. 
We normalize our monic polynomials as 

P n (z) = z n + Pl (n, t)z n ~ l + ... + P n (0), (2.2) 

although we do not always display the t dependence of coefficients of z n ~ l . 

An immediate consequence of the orthogonality condition is the three terms recurrence 
relations; 

zP n (z) = P n+1 (z) + a n P n (z) + /3 n P n _i(z), (2.3) 
with the initial conditions 

P (z) = 1, PoP-i(z) = 0. (2.4) 
An easy consequence of the recurrence relation is 

a n (t) = p 1 (n,t)-p 1 (n + l,t), (2.5) 

and a telescopic sum of the above equation (bearing in mind that p 1 (0, t) = 0), leaves 

n-l 

3=0 

First let us discuss the derivatives of a n and (5 n with respect to t, as this yields the simplest 
equations. We keep wo quite general, as long as the moments, 

jU^(t) := / x 1 wq(x)wj(x; t)dx, % — 0, 1, (2.7) 
Jo 

of all orders exist. Taking a derivative of h n with respect to t gives, 

ti n (t) = -Bw (t)P n (t,tf (2.8) 

or 

(log h n )' = -R n , (2.9) 

where 

R n (t) := Bw (t)^^. (2.10) 
Because f3 n = h n /h n -\ , we have the first Toda equation, 
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Let D n (t) be the Hankel determinant given in (II. 3p . It is well-known that 



n— 1 



D n = Y[hi(t). 



i=0 



In view of ( 12. 9p . it is easily seen that 



Also, 



A n—1 

-\ og D n (t) = -J2Rr (2-11) 

3=0 



d f 1 

= It / P n{ X ) P n-l{x)Wo(x)w ] (x\t)dx 



and therefore 



-Bw {t)P n {t,t)P n ^{t,t) + h n ^ ^(M) 

at 



— Pi(n,f) = 5w (t) =: r„,(t). 



This, together with (12. 5ft . implies 



at 



To summarize we have 



Theorem 3 The recurrence coefficients a n , (3 n and the auxiliary quantities R n , r n satisfy 
the coupled Toda equations, 

(3' n =(R n ^-R n )(3 n , (2.12) 

a' n = r n -r n+1 , (2.13) 

and 

p[(n,t) =r n {t). (2.14) 

Note that the above system (I2.12p - (12.13p does not close, unlike the Toda equations under 
the first time flow. 

It will transpire from next section that we will encounter two more auxiliary quantities 
x n and y n . The ultimate aim of the next section are the expressions of a n and (3 n in terms 
of the auxiliary quantities, and certain identities involving R n , r n , x n and y n . 
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3 Ladder operators, compatibility conditions and dif- 
ference equations 

In this section, we give an account for a recursive algorithm for the determination of the 
recurrence coefficients a n and (3 n , based on a pair of ladder operators and the associated 
supplementary conditions; see (Si), (S 2 ) and (S' 2 ) below. Such operators have been derived 
by various authors over many years. For a quick guide to the relevant literature, the reader 
is directed to the references of [TO] and [6j. In fact Magnus [20] traced this back to Laguerre. 
The derivation where the weight has discontinuities can be found in [I]. 

We find the form of the ladder operators and the compatibility conditions set out below 
well suited for our purpose. 

Let w(x) — wo(x)wj(x; t) . For a sufficiently well-behaved "reference" weight wq(x) (see 
[9] for a precise statement) of the form 

w ( x ) = e -^ x \ 

the lowering and raising operators are 

P' n {z) = -B n {z)P n {z) + p n A n (z)P n ^(z), (3.1) 

P'^z) = \B n {z) + v' {z)]P n - X {z) - A n ^ 1 (z)P n (z), (3.2) 

with 

An{z) :=^ + lf 1 m^M Pl(x)w{y)iv , (3.3) 

Z — I fl n Jo z ~ V 

B n (z) := ^ + ^ T <{z) ~ <{y) P n (y)Pn-i(y)w(y)dy, (3.4) 
z-t h n -\ Jo z-y 

where we have assumed that w (0) = w (l) = 0. Additional terms would have to be included 
in the definitions of A n (z) and B n (z) , if ^o(O), u>o(l) ^ 0. Due to the discontinuity of u>j at 
t, A n (z) and B n (z) have apparent poles at t with residues R n (t) and r n (t) , respectively. 

A direct calculation produces two fundamental supplementary (compatibility) conditions 
valid for all z\ 

B n+ i(z) + B n (z) = {z- a n )A n (z) - v' (z) (Si) 

l + (z- a n ){B n+1 (z) - B n (z)) = (3 n+1 A n+1 (z) - f3 n A n ^{z). (S 2 ) 

We note here that (Si) and (S 2 ) have been applied to random matrix theory in [31]. It 
turns out that there is an equation which gives better insight into the a n and f3 n , if (Si) 
and (S 2 ) are suitably combined. This equation was derived in [10] for smooth weights and 
also holds at the presence of discontinuities. Because its derivation is very simple, we recall 
here the derivation. 

Multiplying (S 2 ) by A n (z) , we see that the right hand side of the resulting equa- 
tion is a first order difference, while the left hand side, with (z — a n )A n replaced by 
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B n+1 (z) + B n (z) + y' Q {z), is a first order difference plus A n (z) . Taking a telescope sum 
together with the appropriate "initial condition" 

B Q (z) = A^{z) = 0, 

produces, 

n-1 

Bl(z) +y' (z)B n (z) + Y,M Z ) = PnA n (z)A n ^(z). (S' 2 ). 

This last equation will be highly useful in what follows. The equations (Si), (S2) and (S' 2 ) 
were also stated in [19J. 

If wq is modified by the multiplication of "singular" factors such as \x — t\ a , then the 
ladder operator relations, (Si), (S 2 ) and (S 2 ) remain valid; see [8]. 

Let ^(z) = P n (z) , eliminating P n _i(z) from the raising and lowering operators gives 

9»(z) - fa) + ^|y) n*) + (b'M - + ^ MZ) ) ^ = °' (3 - 5) 

where we have used (S' 2 ) to simplify the coefficient of *f?(z) in (13.51) . 

For the problem at hand, Wq(x) = x a (l — x)^, we must suppose that a > and (3 > 
so that our weight is suitably well-behaved. Therefore, 

w (x) := x a (l - xf, are [0,1], (3.6) 
M {z) := -a log z - /31og(l - z), 

a P 



v o(^) 



z z — 1 ' 



z-y zy (z-l)(y-l)' 

Substituting (3.7) into the definitions of A n (z) and B n {z) , followed by integration by parts 
gives, 

R n (t) x n {t) x n (t) - R n (t) 

A n {z) = + , 

z — t z — 1 z 

JD ( \ r n(t) _ Vn(t) y n (t) ~ r n (t) - n 

z — t z — 1 z 



where 



ftn Jo 1 -y 

yn{t) := ^_ /' g^zlM _ t)% . 

"•n-l Jo 1 y 
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Note that, 



A n (z) = O 



n ( 1 

B n (z) = — + 0' 



z \z 2 



as z — > oo , which is to be expected. 

We now begin with the characterization of the auxiliary quantities R n , r n , x n and y n . 

Substituting the expressions for A n (z) and B n (z) into (Si) and (S' 2 ) , which are iden- 
tities in z , and equating the residues of the poles at z — 0, 1 and z = t , we find 7 distinct 
difference equations and 3 formulas which importantly performs the sums ^j=o Sj=o 
and Y^ZliRj ■ 

From (S'i), we find, 

r n+1 + r n = (t - «„)#„, (3.8) 

- (yn+l + J/n) = («n ~ l)&n + P, (3-9) 

y n+ i + y n - r n+ i - r n = 2n + 1 + a - a n (x n - R n ). (3.10) 
From (S'2), we find, 

P n R n R n _ x = r 2 w (3.11) 

f3 n x n x n _ x = yl + (3y n , (3.12) 

p n (x n - Rn)(x n -i - = (y n - r n - nf - a(y n -r n - n), (3.13) 

(3 n (x n R n -i + x n -iR n ) = (2n + a + (3)y n - (2n + a)r n + 2y n r n - n(n + a), (3.14) 



E 1 (2n + a + (3)y n - n(n + a) 
3 = 1 

(2n + a + P)(y n - r n ) - n(n + a) 
+ — t , (3-15) 



E" 1 _ (2n + a + /3)(t/» - r») - w(n + a) 
X 3 ~ l-t 

3=0 

+ (2n + a + (3)r n + n(n + a + (3), (3.16) 



- (2n + a + /3)r„ + n(n + a + /5). (3.17) 
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Performing similar calculations on (5*2) , we find, 

(t - a n )(r n+1 - r n ) = (3 n+1 R n+1 - fi n R n - X , (3.18) 

(1 - a n )(y n - y n+ i) = f3 n x n -i - (3 n+ iX n+1 , (3.19) 
- a n (y n+1 - y n + r n - r n+1 - 1) = f3 n+1 x n+1 - (3 n x n -i + (3 n Rn-i - (3 n+1 R n+ i. (3.20) 
Remark 1. It t = , it is clear that r n (0) = R n (0) = 0. So it follows from (I3.14p that 

n(n + a) 



Vn(0) 



2n + a + (3' 



Remark 2. The equation (I3.17P is a consistency check on (13.151) and f)3.16p . 
Remark 3. In [1] and [0], only (Si) and (S' 2 ) are required to provide a complete description 
of the recurrence coefficient, while in [7] and [12], it was found that (S'2) also plays a crucial 
role. We shall see later it is also the case here. 

We now manipulate the equations (I3.8p - (l3.20p with the aim of expressing the recurrence 
coefficients a n and f3 n in terms of r n , R n , x n , y n and of course n and t . We start with 
an easy lemma. 

Lemma 1 We have the equations, 

x n = 2n + l + a + (3 + tR n , (3.21) 
a n = y n+1 -y n + t(r n -r n+1 ), (3.22) 
Pl (n,t) = -y n (t)+tr n (t). (3.23) 

Proof. The equation (13.211) is obtained by summing fl3.8p - fl3.10p . If we substitute the right 
hand side of fl3"TT8"|) and f!3TT9l into fl3^0D . we find <KT3 . The equation f!3T2^D follows from 
a telescopic sum of f)3.22p . keeping in mind the "initial" conditions, r (t) = yo(t) = 0, and 
EjTo a 3 = ~Pi(n, t) ; see (EU). □ 

Differentiating 03.231) with respect to t and taking into account (12.141) . we find in the 
next lemma 

Lemma 2 

Vn = K- (3-24) 

Now, we are ready to express the recurrence coefficients a n and (3 n in terms of auxiliary 
quantities. In particular, we emphasize here the importance of p x (n, t) as stated in the 
following theorem. 

Theorem 4 The recurrence coefficients a n and (3 n have the following representations: 

(2n + 2 + a + (3)a n = 2tr n - 2y n - (3 + (2n + 1 + a + (3)t + (1 - t)x n , 

= 2 Pl (n, t)-P+(2n+l + a + P)t+(l- t)x n . (3.25) 
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and 

(2n - 1 + a + (3){2n + 1 + a + /3)/3 n 

= (y n - tr n ) 2 + (2n + a)tr n + [f3 - (2n + a + /?)t]y n + n(n + a)t 
= [p^n, t)] 2 + (2n + a)p 1 (n, t) + (2n + a + /3)(1 - t)y n + n(n + a)t 
= [ Pl (n, t)] 2 + + (2n + a + (3)t] Pl {n, t) 

+ (2n + a + f3)t(l - t)p[(n,t) +n{n + a)t. (3.26) 

Proof: First we eliminate r n+ \ from (12.131) and (13. 8p to find, 

2r n = (t - a n )i^ + <. (3.27) 

In view of ( 12.131) . we can replace r n — r n+ i in ( 13. 22ft by o4 , i.e., 

a n = Vn+i -Vn + ta' n . 

Now eliminate y n+ \ from above and (13.91) . we obtain 

2y n = -a n + ta' n - (3 + (1 - a n )x n . (3.28) 

The first equality in f!3.25j) then follows from further eliminating a' n (t) by (I3.27P and (13.281) 
and the second equality is immediate on account of (13.231) . 

To derive (13.261) . we use (13.211) to eliminate x n and x n _\ in (13.121) . and make use of 
(13~TTD to find, 

Pn\pnCn—l tC n R n —i -\- tC n —\R n -\- t R n _\R n \ 
— Pn[ c n c n-1 + t c nRn-l + tC n _iR n ] + tV 2 

= y 2 n + f3y n . (3.29) 

where c n := 2n + 1 + a + (3 . 

On the other hand, it follows from (EDI) , (EOT!) and (133TD that 

Pni^nRn—l •^n—lRn) 

n-1 + ( c n-l + tR n -i)R n ) 

= (3 n {c n R n _i + c n _ii4) + 2tr 2 

= (2n + a + /3)y n - (2n + a)r n + 2y n r n - n(n + a), (3.30) 

or equivalently, 

ftniC-nRn—l Cn—lRn) 

= (2n + a + P)y n - (2n + a)r n + 2y n r n - 2tr 2 n - n(n + a). (3.31) 

Substitution of the above formula into ( I3.29P gives us the first equality of (13.261) . The last 
two equalities follow from (I3.23P and (12.141) . □ 
Recall that H n (t) is the logarithmic derivative of D n (t) defined in (II. 5p . The next 
theorem contains an important identity expressing H n (t) in terms of the auxiliary quantities 
y n and r n . This will in turn give us expressions of y n and r n in terms of H n (t) and H' n (t). 
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Theorem 5 The logarithmic derivative of the Hankel determinant (11.51) can be expressed in 
terms of y n and r n as 



H n (t) = (2n + a + (3)(y n - tr n ) - n(n + a) 
= — {2n + a + /3)p 1 (n, t) — n{n + a), 



while 



r n (t) 

Vn(t) 



2n + a + (3 
-tH' n {t) + H n {t) + n{n + a) 
2n + a + (3 



(3.32) 

(3.33) 
(3.34) 



Proof: Combining (12 .1 If) . which gives the logarithmic derivative of D n (t) as the sum 
5^ ■ Rj , and (13.151) . which gives the sum in terms of y n and r n , leaves, 



d_ 

It 



In IX 



(2n + a + (3)(y n — r n ) — n{n + a) {2n + a + (3)y n — n(n + a) 



1-t 



After a little simplification, equation (13.321) follows. Taking a derivative with respect to t 
on the second equality in (13.321) . it follows from (12. 14j) that 

H' n (t) = -(2n + a + P)r n {t), 

which is (13.331) . The derivation of (13.341) is now obvious. 
Remark 4. With reference to (I3.23|) and Remark 1, we see at t = 0, 

n(n + a) 



□ 



p x (n, 0) 



2n + a + f3' 



Integrating (3.33), we find, 



D n (t) = D n (0) exp 



{2n + a + (3) [ 
Jo 



Pi(n, s) - pj(n,0) 



ds 



s(l -a) 

where D n (0) is Hankel determinant associated with "unpertured" Jacobi weight 

{A + B)x a (l-xf, 

and it is explicitly known. 

Following what was done in [3], [7], [12] and [6], we derive two alternative representations 
of R n in terms of y n and r n , which is stated in the next proposition. 
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Proposition 1 The auxiliary quantity R n has the following representations 

R j f) = (2n + l+a + (3) [l(r n , y ni t) - (1 - t)y' n (t)] 

2k(r n ,y n ,t) 



1 l(r n ,y n ,t) + {l-t)y' n {t) 



(3.36) 



R n (t) 2{2n + 1 + a + /3)r* ' 
where 

l(r n , y n , t) := (2n + a + (3)y n - (2n + a)r n - 2tr 2 n + 2y n r n - n(n + a) (3.37) 

and 

k(r nj y n ,t) := (y n - tr n ) 2 + (2n + a)tr n + {(3 - (2n + a + f3)t)y n 

+ n(n + a)t (3.38) 

Proof: Using (13. lip , we eliminate R n -\ in (I3.3ip and obtain 

2 

(2n + 1 + a + p)^- + (2n - 1 + a + (3)(3 n R n 

Rn 

= (2n + a + 0)y n - (2n + a)r n + 2y n r n - 2tr 2 n - n(n + a). 

Replacing (3 n in the above formula with the aid of (13.261) . we have 

(2n + 1 + a + (3)r 2 k(r n ,y n ,t) 

R n 2n + l + a + p n 

= (2n + a + {3 + 2r n )y n - (2n + a)r n - 2tr 2 n - n{n + a). (3.39) 

This is a linear equation in R n and 1/R n . We shall derive another linear equation in R n 
and 1/Rn- 

To proceed we first note the first Toda equation (I2.12p . reproduced here; 

fln = (Rn-l — Rn)(3 n 

Applying | to the first equality in fl3T26l) . it follows from flSTTTD . (13T26D . fl3T2ip and the above 
formula that 

,2 

(2n - 1 + a + (3) (2n + 1 + a + /3)^- - fc(r n , y ft , t)i? n 

= 2tr 2 n + (2n + a- 2y n )r n - (2n + a + /3)y n 
+ (2n + a + /5)(l-t)y; + n(n + a). (3.40) 

The equations (13.351) and (13.361) now follow from solving for R n and l/i?„ from (I3.39P and 
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4 Proof of Theorem H and derivation of Painleve VI 

With the above preparations, we are now ready to prove Theorem [1] and show that x n (t) is, 
up to a linear fractional transformation, a particular Painleve VI. 
Proof of Theorem [I]: Multiplying (I3.35j) and (I3.36j) gives us 

(1 - t) 2 [y' n (t)] 2 = l 2 (r n , y n , t) - Ak(r n , y n , t)r 2 n , (4.1) 

where l(r n ,y n ,t) and k(r n ,y n ,t) are given by (13.371) and f!3.38j) . respectively. Substituting 
(I3.33P and (I3.34j) into (14.11) yields a second order non-linear differential equation for H n . 
Using (11.61) to replace H n by a{t) , we finally get (ll.7p . 

The initial conditions follow from (13.331) and the fact that r n (0) = 0. This completes 
the proof of our theorem. □ 

The rest of this section is devoted to the derivation of a differential equation satisfied 
by x n . As this is going to quite complicated algebraically, we will separate the proof into 
several parts. 

Proposition 2 The auxiliary quantity x n has the following representations 

x ,q = (2n+l + a + {3) [7(r n , y n , t) - t(l - t)y' n (t)} ^ 

2k(r n , y n , t) 
1 T(r n ,y n ,t)+t(l-t)y' n (t) 



x n (t) 2{2n + 1 + a + (3)((3 + y n )y n 
where 

T(r n , y n ,t) := 2y 2 n + {2(3 - (2n + a + (3 + 2r n )t)y n + (2n + a)tr n + n(n + a)t (4.4) 
and k(r n ,y n ,t) is given in (13.3 



Proof: The equality (14.21) follows from a combination of (13.211) and (13.351) . To show (14.31) . 
we use (I3.2ip to eliminate R n and R n -i in (13.111) and get 

Pn^n^n—l Pn^-n—l'^n f^rfirfin—\ t ^ ni (^"^) 

where c n := 2n + 1 + a + /3 . On account of (13.121) . (13.261) and (I3.38j) . the above formulas is 
equivalent to 

2 , n c n(Vn + fiVn) k(r n ,y n ,t)x n . 2 2 t a r \ 

Vn + PVn + k{r n , y n , t) = t r n . (4.6) 

Xn Cn 

Solving gS]) with the aid of fl4~2]) yields (Ojl . □ 
Next, we express r n in terms of x n , x' n and y n . 
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Differentiating the first formula in f!3.25j) with respect to t , it follows from (I2.13P and 
(153311 that 

(2n + 2 + a + /3)(r n - r n+1 ) = 2r n - x n + (1 - t)< + 2n + 1 + a + (3. (4.7) 

By ( 13.81) and (13.211) . it is also easily seen that 

(t - a n )(x n - 2n - 1 - a - (3) = t(r n + r n+1 ). (4.8) 

Eliminating r„ +1 from the above two formulas, together with the first equality of (13.251) . we 
obtain an expression of r n in terms of x n , x' n and y n , in the next Lemma, 

Lemma 3 

I 2n + 1 + a + (3 + (l-t)x' n (2y n + (3-(l-t)x n + t)(2n + l + a + (3-x n ) 



2 2^c^ 2/>^c yi 

(4.9) 

As the final step involves rather large formulas, we will only give a qualitative discussion 
here. 

We begin by substituting (14.91) into (14. 2p and (14. 3p . After some computations, we find a 
pair of linear equations in y n and y' n , with coefficients depending t . Solving; 

this linear system gives us 

Vn{t) F(x n , X n , 71, t), 

y'n(t) = G(x n ,x' nl n,t), 

where the form of F and G are explicitly known. We note here x" n does not appear either 
in F or in G . Since the expressions are somewhat long, we are not prepared to write them 
down. 

From the obvious equality 

-^F{x ni x n} n } i) = y n {t) = G{x n} x n} n } t') } 

it can be shown, after some computations, that the above equation is equivalent to 

{{2n + a + (3){2n + a + (3 + 1) + \{2n + a + (3 + l)t - 2(2n + a + (3) - l]x n (t) 

- (f - l)x 2 n (t) + t(t - l)x' n (t)}$(x n , x' n , xl n, t) = 0, (4.10) 

where $ is a functions that is explicitly known. Clearly, the above formula yields two 
differential equations. One is a Riccati equation, whose solutions are of no relevance to our 
problem and should be discarded. 
If we make the substitution 

x n (t) = (2n + l + a + P) l ~^_ n ® (4.11) 
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in $ = , we find after some computations that, W n (t) satisfies the following Painleve VI 
equation 



2\W n W n -1 W n -tJ y n/ \t t-1 W n -t, 
Wn (W n -l)(W n -t) f bt_ c(t-l) dt(t-l) \ 

t 2 (t - 1) 2 \ a+ wy (w n -iy + (w n -ty) l4 - iZj 

with 

(2n + a + [3 + l) 2 , a 2 (3 2 , 1 

a = , o = , c = — , a = - 

2 ' 2 ' 2 ' 2 

and initial conditions 

W n (0) = 0, <(0) = 1, 

as advertised in the introduction. The above initial conditions follow from (13.211) and the 
fact R n (0) = 0. 



5 Proof of Theorem [2 



We begin with the integral representation of D n (t) given in fll.4f) . which in the situation 
where A = and 5 = 1 reads 

If n 

D n (t) = - [A n (x)] 2 T\x%(l-x k ) p dx k , A n (x)= TT (xj-Xi). (5.13) 

"MM] k=l l<«<i<n 

A simple change of variable Xk = Uk + t in the above multiple integral yields 

If n 
D n (t) = - / [A n (y)] 2 T]>* + - t - y*)%* 
n! ^[0,l-i]» t=i 

= (1 _ t) n(n + P) / [An(y)] 2 JT [(1 _ t) + ^ (1 _ ^ 
™ ! i[o,i]« 

where z — 1 — As £ — > 1~, we find, 

1 /" n 
A,(f) ~ (1 - t) n ^ - / [A n (j/)] 2 TT(1 - y k fdy k . (5.15) 

To evaluate the multiple integral on the right hand side of (15.151) . we set 

If n 
D n [w a , {-1, 1)] := - / [A n (x)} 2 TT(1 + x k ) a (l - x k fdx k . (5.16) 
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Then, an elementary translation and scaling shows that 

1 f n 
D n K^(-l,l)] = 2< n+ « + V - / [A n (y)] 2 n^(l -Vkfdy k 

= 2 n( - n+a+ ®D n [w a>p {0,l)) (5.17) 

By [3, (1.6)], we have 

D n [w a , p {-\, 1)] = 2-< n+a+ ^ (2vr)" F(a, (3) K(a, (3, n), (5.18) 

where F(a,/3) and K(a,j3,n) are given in (j 1.9 ft and fll.lOp . respectively. Hence, it follows 
from the above two formulas that 



D n [w a A0A)} = I [A n (y)] 2 TT^(l -y k fdy k 

= 2- 2 < n+a+ ^ (27r) n F(a, P)K(a, (3, n). (5.19) 
Note that D n (0) = D n [w a ^(0, 1)] , combining floTTD and fl5TT9"|) gives 

D n(t) Dn[w O ,p(0,l)] (1 _ t) n(n+f3) 

D n (0) ~ D n [w a , (O,l)V 

F(a, p)K(a, p, n) 

which is fOjl . 

This finishes the proof of Theorem [2j □ 
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